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Superconducting Puddles and "Colossal" Effects 
in Underdoped Cuprates 
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Phenomenological models for the antiferromagnetic (AF) vs. d-wave superconductivity competi- 
tion in cuprates are studied using conventional Monte Carlo techniques. The analysis suggests that 
cuprates may show a variety of different behaviors in the very underdoped regime: local coexistence 
or first-order transitions among the competing orders, stripes, or glassy states with nanoscale su- 
perconducting (SC) puddles. The transition from AF to SC does not seem universal. In particular, 
the glassy state leads to the possibility of "colossal" effects in some cuprates, analog of those in 
manganites. Under suitable conditions, non-superconducting Cu-oxides could rapidly become su- 
perconducting by the influence of weak perturbations that align the randomly oriented phases of 
the SC puddles in the mixed state. Consequences of these ideas for thin-film and photoemission 
experiments are discussed. 

PACS numbers: 74.20.De,74.72.-h,74.20.Rp 



I. INTRODUCTION 



Clarifying the physics of high-temperature supercon- 
ductors (HTS) is still one of the most important chal- 
lenges in condensed-matter physics. There is overwhelm- 
ing experimental evidence for several unconventional 
regimes in HTS, including a pseudogap region at tem- 
peratures above the superconducting (SC) phase, and 
a largely unexplored glassy state separating the parent 
antifcrromagnet (AF) from the SC phase at low hole- 
doping x. In addition, recent investigations unveiled an- 
other remarkable property of HTS's that defies conven- 
tional wisdom: the existence of giant proximity effects 
(GPE) in some cuprates,[H H Q where a supercurrent 
in Josephson junctions was found to run through non- 
SC Cu-oxide-based thick barriers. This contradicts the 
expected exponential suppression of supercurrents with 
barrier thickness beyond the short coherence length of 
Cu-oxides. The purpose of this paper is to propose an 
explanation based on a description of the glassy state as 
containing SC puddles. This nanoscale inhomogeneous 
state leads to colossal effects in cuprates, in analogy with 
manganites 0.0.0 In addition, it is argued that differ- 
ent inhomogeneous states could be stabilized in different 
Cu-oxides, depending on coupling and quenched disor- 
der strengths. In fact, neutron scattering studies have 
revealed "stripes" of charge in Nd-LSCO,0,@ but scan- 
ning tunneling microscopy (STM) experiments^. ITo| in- 
dicate "patches" in Bi2212, consistent with our analysis. 
There is no unique way to transition from AF to SC. 

Studies of the t-J model have revealed SC and striped 
states^], ^3 evolving from the undoped limit. Then, it 
is reasonable to assume that AF, SC, and striped states 
are dominant in cuprates, and their competition regu- 
lates the HTS phenomenology. However, further compu- 
tational progress using basic models is limited by clus- 
ter sizes that cannot handle the nanoscale structure un- 



veiled by STM experiments. Considering these restric- 
tions, here a phenomenological approach will be pursued 
to understand how these phases compete, incorporating 
the quenched disorder inevitably introduced by chemical 
doping. This effort unveils novel effects of experimen- 
tal relevance, not captured with first-principles studies. 
Two models are used, one with itinerant fermions and 
the other without, and the conclusions are similar in 
both. Hopefully, this effort will jump start a more de- 
tailed computational analysis of phenomenological mod- 
els in the high-Tc arena, since most basic first-principles 
approaches, including Hubbard and t — J investigations, 
have basically reached their limits, particularly regarding 
cluster sizes that can be studied. 



II. MODEL I: ITINERANT FERMIONS 

The analysis starts with a phenomenological model of 
itinerant electrons (simulating carriers) on a square lat- 
tice, locally coupled to classical order parameters: 

H F = -t (cL c j^ + ^- c -) + 2^Ji^fsf-^Min iCT 

<ij>,fT i icr 



where c- la are fermionic operators, sf=(n^ — n^)/2, n- la 
is the number operator, D is the lattice dimension, and 



Air 



are complex numbers for the SC order 



parameter defined at the links (i,i+a) (a = unit vector 
along the x or y directions). At Ji=0, d-wave SC is fa- 
vored close to half-filling since the pairing term involves 
nearest-neighbor sites, as in any standard mean-field ap- 
proximation to SC. The spin degrees of freedom (d.o.f.) 
are assumed to be Ising spins (denoted by S?). Stud- 
ies with 0(3) d.o.f. were found to lead to qualitatively 
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FIG. 1: (a) MC phase diagram for Eq. Q without disorder at 
low temperatures. Instead of presenting a three dimensional 
phase diagram we have chosen to present a two dimensional 
cut along V=l-J/2 for simplicity. Five regions are observed: 
AF, d-SC, stripes, coexisting SC+AF, coexisting stripes+SC, 
and metallic, (b) MC phase diagram including temperature 
along "Path 1". (c) MC phase diagram along "Path 2". Lat- 
tice sizes in all cases are 8x8 and 12x12. (d) n vs. n along 
Paths 1 and 2. Transitions along Path 1 appear continu- 
ous, whereas along Path 2 there are indications of first-order 
transitions, (e) Spin structure factor S(q) at (tt,tt) and for 
incommensurate (IC) momenta. 



similar conclusions, but they are more CPU time con- 
suming. The parameters of relevance are J\, /Zi, and 
Vi (t is the energy unit), that carry a site dependence 
to easily include quenched disorder which is inevitable 
in chemically doped compounds as the cuprates. For a 
fixed configuration, {Aj Q } and {S?}, the one-particle sec- 
tor is Bogoliubov diagonalized. In the limit T^0, the 



Bogoliubov-de Gennes equations are recovered minimiz- 
ing the energy 00 0. Then, a standard Monte Carlo 
(MC) simulation similar to those for Kondo-lattice mod- 
els is carried out (details in Ref. |4|). One of the goals is 
to estimate T c , as well as T*, roughly defined as the tem- 
perature at which strong short-distance SC correlations 
develop (more details are given below) . Finally, note that 
model Eq. JIJ is not derived but proposed as a possible 
phenomenological model for AF vs. SC competition. The 
results will be shown to justify this assumption. More- 
over, the qualitative simplicity of our conclusions sug- 
gests that similar models will lead to similar scenarios. 



A. Phase Diagram in the Clean Limit 

Without quenched disorder, Vi, J\ and /ii in Eq. (JJJ are 
site independent. The standard MC analysis carried out 
in these investigations (details provided below) reveals 
that in the clean limit the low temperature (T) phase 
diagram, Fig^a), has a robust AF phase for electronic 
densities n~l and a d-wave SC phase for n<l. The d- 
wave correlation function, defined as 



C s f (m)=£(|A 1 ||A i+m |co S (C- 



1 )}, ( 2 ) 



was used to estimate T c as the temperature at which 
d-wave correlations at the largest distances for the lat- 
tices considered here are 5% of their maximum value at 
|m|=0. The 5% criterion is arbitrary but other crite- 
ria lead to identical qualitative trends, slightly shifting 
the phase diagrams. T* is deduced similarly, but us- 
ing the shortest non-zero distance correlation function 
(|m|=l). The Neel temperature, In, associated with the 
classical spins was defined by the drastic reduction (< 
5% of | m |=0 value) of the long-distance spin order using 
Cs(m)=2j (5f5f +m ), while relates to short-range 
spin order. The results presented in Fig. ^a) are not sur- 
prising since these states are favored explicitly in Eq. 
by the second and fifth terms, respectively. However, 
the phase diagram presents several nontrivial interesting 
regions: (i) Along "Path 1" in Fig.Ufa), the AF-SC tran- 
sition occurs through local coexistence, with tetracritical 
behavior (Fig.ID»).0 (ii) Along "Path 2" the AF-SC 
interpolating regime has alternating doped and undoped 
stripes (stripes in MC data are deduced from spin and 
charge structure factors, and low-T MC snapshots), and 
a complex phase diagram Fig.^c). These stripes evolve 
continuously from the V=0 limit that was studied be- 
fore by Moreo et al., and as a consequence we refer the 
readers to Rcf. for further details on how stripes were 
identified. It remains to be investigated if these stripes, 
involving SC and AF quasi-lD lines, have the same or dif- 
ferent origin as those widely discussed before in the high- 
Tc literature.000|2al2l|22 At V ^0, the doped 
regions of the stripes have nonzero SC amplitude at the 
mean- field level. |23| In view of the dramatically different 
behavior along Paths 1 and 2, we conclude that in our 



model there is no unique AF^SC path. This is in agree- 
ment with experiments since La2_ x Sr a; Cu04 (LSCO) and 
others have stripes, 0, H, |24| while Ca^NayCuO^Cl? 
has a more complex inhomogeneous pattern. 10] Both, 
however, become SC with increasing x. This suggests 
that the underdosed region of Cu-oxides may not be uni- 
versal. 



B. Phase Diagram with Quenched Disorder 

Our results become even more interesting upon in- 
troducing quenched disorder, with a MC phase diagram 
shown in Fig. b) . The similarity with the widely ac- 
cepted phase diagram of the cuprates is clear. The dis- 
order has opened a hole-density "window" where none 
of the two competing orders dominates. Inspecting "by 
eye" the dominant MC configurations (snapshots) at low- 
T in this intermediate regime reveals a patchy system 
with slowly evolving islands of SC or AF, and random 
orientations of the local order parameters, leading to an 
overall disordered "clustered" state. In Fig.Etb), a new 
temperature scale T* at which the fermionic density-of- 
states (DOS) develops a pseudogap (PG) (Fig. Et c )) was 
also unveiled. The AF and ci-SC regions both have a 
gap (smeared by T and disorder, but nevertheless with 
recognizable features). But even the "disorder" regime 
(case b in Fig. Etc)) has a PG. MC snapshots explains 
this behavior: in the disordered state there are small 
SC or AF regions, as explained above. Locally each has 
a smeared-gap DOS, either AF or SC. Not surprisingly, 
the mixture presents a PG. The behavior of T* vs. x 
is remarkably similar to that found experimentally. The 
cuprates' PG may arise from an overall-disordered clus- 
tered state with local AF or SC tendencies, without the 
need to invoke other exotic states. This PG is corre- 
lated with robust short-range correlations (dashed lines 
in FigEtb), see caption for details.). 

The numerical procedure that led to Figs. ^ an d El ls 
standard, well-known in the manganite context where 
formally similar models are widely studied, 0, 0, thus 
here only a few representative results will be presented 
for completeness. For instance, Fig.[3ta) shows the order 
parameters along "Path 1" of Fig. ^a), indicating how 
each ordered region was determined in the clean limit. 
Clearly, a region of coexistence can be observed. Similar 
data were used to complete the phase diagram. Likewise, 
Fig. Etb) contains the order parameter vs. T with and 
without disorder, in the interesting region of couplings 
and doping. There is a drastic difference between clean 
and dirty limits, the latter showing no global dominant 
order. However, examining relevant MC configurations, 
small SC and AF clusters with random orientations are 
found. 
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FIG. 2: (a) Plaquette impurity schematic representation. 
Disorder may have several forms, but here we mimic Sr-doping 
in single-layers. Sr 2+ replaces La 3+ , above the center of a Cu- 
plaquette in the Cu-oxide square lattice. Then, as hole carri- 
ers are added, a hole-attractive plaquette-centered potential 
should also be incorporated. Near the center of this potential, 
n should be sufficiently reduced from 1 that, phenomenolog- 
ically, tendencies to SC should be expected. To interpolate 
between the SC central plaquette and the AF background, a 
plaquette 'halo' with no dominant tendency was introduced 
Parameters are chosen such that the blue (black) region favors 
superconductivity, (J, V, p) = (0.1, 1.0, —1.0), with a surround- 
ing white region where (J, V, /x)=(0.1, 0.1, —0.5) with no order 
prevailing. The impurity is embedded in a background (red, 
dark gray) that favors the AF state, (J, V, /i)=(1.0, 0.1, 0.0). 
However, the overall conclusions found here are simple, and 
independent of the disorder details, (b) MC phase diagram 
for model Eq. Q including quenched disorder (lattices stud- 
ies are 8x8 (results shown) and 12 x 12). Shown are T c and Tn 
vs. number of impurities (equal to number of holes). The SC 
and AF regions with short-range order (dashed lines), and T* 
as obtained from the PG (dot-dashed line) are also indicated, 
(c) DOS at points a, b, and c of (a), with a PG. 



III. MODEL II: LANDAU GINZBURG 

The results reported thus far, based on Eq. (Q, have 
already revealed interesting information, namely the pos- 
sible paths from AF to SC, and a proposed explanation of 
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FIG. 3: (a) MC S(n,n) and superconducting correlation at 
maximum distance, Cf c y , Eq. along "Path 1" (Fig. 1(a)), 
using an 8 x 8 lattice and at low temperature T=0.02. (b) AF 
and SC correlations at maximum distance for the model with 
8 impurities and without disorder (clean). The latter is the 
point J=0.6 of "Path 1". The disordered case corresponds 
to 8 plaquettes on a 12x12 lattice. Typically, 5,000 sweeps 
were used for thermalization and for measurements. With 
quenched disorder, many points were obtained after averaging 
over JVdis^lO disorder realizations. The results were found 
to mildly depend on the disorder configuration, thus many 
results were obtained with a smaller iVdis- 



the glassy state as arising from the inevitable quenched 
disorder in the samples. However, the inhomogeneous na- 
ture of the clustered region suggests that percolative phe- 
nomena may be at work, and larger clusters are needed. 
To handle this issue, another model containing only clas- 
sical d.o.f. is proposed, with low-powers interactions typ- 
ical of Landau- Ginzburg (LG) approaches: 



+a 



H = n J2 lAil 2 + f £ |Ai| 4 + 5>(i, a)S, ■ S, 

i i i.a 

- Pi(h a)|Ai||A 1+0 | cos(*i - *;+«) + r 2 £ |Si| 2 

i,a i 

+ ^j^r+^Ei^iw- 0) 



The fields Aj=| Aj|e ' are complex numbers repre- 
senting the SC order parameter. The classical spin 
at site i is Si=|Si|(sm(0i)cos(0i),sm(0i)sin(<^i),cos(0i)). 
pi(i, a) = l — p 2 (i, a) is used as the analog of V—l-J/2 of 
the previous model to reduce the multiparameter charac- 
ter of the investigation, allowing an AF-SC interpolation 
changing just one parameter, a denotes the two direc- 
tions x and y in 2d, and also z for multilayers. p"i$, ct) 
was chosen to be isotropic, i.e., a- independent. 



A. Basic Properties 

Clearly, the lowest-energy state for p2=0 is a homo- 
geneous SC state (if a)—pi>0). When pi=0 the 
lowest-energy state is AF (if p 2 (i, a)—p2>0). In the clean 



limit, this model was already studied in the SO (5) con- 
text, where the reader is referred for further details. Our 
approach without disorder has similarities with SO (5) 
ideas [13 where the AF/SC competition as the cause of 
the high-T c phase diagram was extensively discussed al- 
though nowhere in our investigations we need to invoke a 
higher symmetry group. The relevance of tetracriticality 
in LapCuO^+a has also been remarked by E. Demler et 
q^.pfij and Y. Sidis. 26] In the present work, disorder is 
introduced by adding a randomly selected bimodal con- 
tribution, i.e. p2(i, cv)=p2±VF, where W is the disorder 
strength (W=0 is the clean limit). It is expected that 
other forms of disorder will lead to similar results. 



B. Phase Diagram 

Monte Carlo results for Eq. are in Fig. QJi, for 
"weak" coupling tti2=0.7, which leads to tctracritical be- 
havior. Both at W=0 and W^0, the qualitative simi- 
larity with fermionic model results (Fig's^b) andE^b)) 
is clear. Coexisting SC and AF clusters appear in MC 
snapshots (not shown) . Then, both models share a simi- 
lar phenomenology, and Eq. can be studied on larger 
lattices. The only important difference between the two 
models is that Eq. © cannot lead to doped-undoped 
stripes, but the more general case Eq. does. Fig.0Jb) 
illustrates how the phase diagram, Fig. Ufa), was ob- 
tained. For completeness, note that increasing the cou- 
pling U12 a first-order SC-AF transition can be obtained. 
However, the addition of disorder leads to a very similar 
phase diagram as in the case of U12 = 0.7. This is shown 
in Fig. 2t and is the equivalent of Fig. [3Ji in the regime 
of "strong" coupling. 

Some of the experimental predictions related with our 
SC-AF clustered state are simple (the most elaborated 
ones are in the next section) . In most ways a very under- 
doped cuprate can be tested, there should be two compo- 
nents in the data. For instance, a typical photoemission 
spectra in our framework should have two clearly distinct 
coexisting signals. This result, which will be discussed in 
more detail in a future publication, is compatible with 
photoemission experiments for ir=0.03 LSCO, that re- 
veal spectral weight in the node direction of the d-wave 
superconductor even in the insulating glassy regime. [53 
Nodal d-wave SC particles surviving to low x was ob- 
served in Ref . [28l 



IV. COLOSSAL EFFECTS IN CUPRATES 

One of the main results of these investigations is that 
the models studied here can present "colossal" effects, 
very similarly in spirit as it occurs in manganites. Con- 
sider a typical clustered state (Fig. EJb)) found by MC 
simulations in the disordered region. This state has pre- 
formed local SC correlations - nanoscale regions hav- 
ing robust SC amplitudes within each region, but no SC 
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FIG. 4: (a) MC phase diagram (for Eq. ©) at 
Mi2=0.7. Parameters are n=-l, rz =— 0.85, and ui=it2=l 
but the conclusions are not dependent on coupling fine- 
tuning. Spin C sp m(m) = -i- £\ (Si • Si +m ) and SC correla- 
tions C S c(m) = i J2i |Ai||A i+m | cos(*i - *i+m) were mea- 
sured. The behavior of these functions at the largest (short- 
est) distance determine T c and Tn (T*) (same criteria as 
for Eq. Q). With disorder, the phase diagram (shown) 
has an intermediate "clustered" state with short-range order. 
T* is also indicated (dashed line). Note the similarity with 
Fig.[5Jb). Inset: results at VK=0 showing tetracriticality (ma- 
genta (dark) indicates SC-AF coexistence), (b) AF and SC 
correlations at maximum distance for the model Eq. Q with- 
out and with disorder (W=0.0 and 0.7, respectively). pi=0.5 
and ui2=0.7 were used, using a 24x24 lattice. Typically, 
for the LG model 25,000 sweeps were used for thermalization 
and measurements, (c) MC phase diagram of model Eq. 
at ui2=2. The clean case (W=0, solid lines) is bicritical-like, 
but with disorder 11^=0.5 a clustered region between SC and 
AF opens as well. 



phase coherence between different regions - rendering the 
state globally non-SC (the averaged correlation at the 
largest distances available, Cg^ x , is nearly vanishing). 
Let us now introduce an artificial SC "external field" , 
which can be imagined as caused by the proximity of a 
layer with robust SC order (e.g., comprised of a higher- 
T c material). In practice, this is achieved in the calcula- 
tions by introducing a term |A§q | J2i Pi(h cos(^i), 




FIG. 5: Left: Q n c ax vs. A§£ (see text) on a 24x24 lattice, 
with wi2=0.7 and W=0.5, at the five points a-f indicated in 
Fig. A "colossal" effect is observed in a and 6 where 

the Ag^=0 state is "clustered". A much milder (linear) ef- 
fect occurs far from the SC phase (e and /). MC snapshots 
are shown at Agc'=0.0 right, top and Agc'=0.2 right, bottom, 
both at T=0.1 and p2=0.5, using the same quenched-disorder 
configuration. The color convention is explained in the circle 
(colors indicate the SC phase, while intensities are propor- 
tional to -Re(Ai)). The AF order parameter is not shown. 
The multiple-color nature of the upper snapshot, reflects a 
SC phase randomly distributed (i.e. an overall non-SC state). 
However, a small external field rapidly aligns those phases, 
leading to a coherent state (bottom). 



where Ag|S acts as an external field for SC. The depen- 
dence of CgQ X with Agp is simply remarkable (Fig.jsja)). 
While at points e and /, located far from the SC re- 
gion in Fig.QJi, the dependence is the expected one for a 
featureless state (linear), the behavior closer to SC and 
small temperatures is highly nonlinear and unexpected. 
For example, at point a, Cc^f* vs. AgjS has a slope (at 
Agp=0.02) which is ~250 times larger than at e (^13 
times larger than at W=0, same T, p2, and 1*12-). 

The reason for this anomalous behavior is the clustered 
nature of the states. This is shown in the state Fig.OJc), 
contrasted with (b), where a relatively small field - in 
the natural units of the model - nevertheless led to a 
quick alignment of SC phases, producing a globally SC 
state, as can be inferred from the uniform color of the 
picture. Having preformed SC puddles vastly increases 
the SC susceptibility. Since Fig. Ufa) was obtained in 
a trilayer geometry it is tempting to speculate that the 
proximity of SC layers to a non-SC but clustered state, 
can naturally lead to a CPE over long distances, as ob- 
served experimentally in a similar geometry. Q, Q 
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FIG. 6: (a) MC phase diagram (for Eq.©) at ui 2 =0.7. Pa- 
rameters are n=-l, T2 =— 0.85, ui=W2=l, W = 0.5 with one 
layer (solid colors) and two layers (dashed line). The addi- 
tion of an extra layer increases the critical temperature of the 
superconductor as well as the Neel temperature, (b) T c vs. 
N e for Ui2=0.7, p 2 =0.3, W=0.7, and 24 2 xiV i clusters. Shown 
are results with and without disorder, (c) The experimental 
T c (in K) is shown for three HTS families, as indicated, up to 
3 layers (data from Ref . l29f) . 



V. DEPENDENCE OF T c WITH THE NUMBER 
OF LAYERS 

The nanoscale clusters also leads to a proposal for ex- 
plaining the rapid increase of T c with the number of 
Cu-oxide layers Nt, found experimentally, at least up to 
3 layers. In this effort, the MC phase diagrams of single-, 
bi-, and tri-layer systems described by Eq. ©, with and 
without disorder, were calculated using exactly the same 
parameters (besides a coupling p2{\,z), equal to those 
along x and y, to connect the layers). It was clearly ob- 
served that the single layer has a substantially lower T c 
than the bilayer. This can be understood in part from 
the obvious critical fluctuations that are stronger in 2D 
than 3D. But even more important, cluster percolation 
at Wj^O is more difficult in 2D than 3D (since otherwise 
2D disconnected clusters may become linked through an 
interpolating cluster in the adjacent layer). Then, in the 
phenomenological approach presented here it is natural 



that T c increases fast with Ne, when changing from 1 to 
2 layers as shown in Fig. This concept is even quan- 
titative - up to a scale - considering the similar shape of 
T c vs. Ne found both in the MC simulation and in exper- 
iments (see FigsEjD-c. Note that the subsequent decrease 
of T c for 4 or more layers observed experimentally could 
be caused by inhomogeneous doping, beyond our model). 
Our MC results suggest that the large variations of T c 's 
known to occur in single-layer cuprates can be attributed 
to the sensitivity of 2D systems to disorder. As Nt in- 
creases (the system becomes more 3D), the influence of 
disorder decreases, both in experiments|30j and simula- 
tions. 



VI. CONCLUSIONS 

Summarizing, here simple phenomenological models 
for phase competition showed that - depending on de- 
tails - different cuprates could have stripes, local coexis- 
tence, first-order transitions, or a glassy clustered state 
interpolating between AF and SC phases. Figure [7| il- 
lustrates our proposed possibilities. In Cu-oxides where 
the glass state is realized, namely where SC puddles are 
present, this study revealed the possibility of colossal ef- 
fects. A schematic representation of the proposed glassy 
state with colossal effects is in Fig. [S] This proposal 
could provide rationalization of recent results in trilayer 
thin-film geometries. Q, 0, 

After submission of this work, we learned of interest- 
ing experimental efforts that complement the discussion 
presented here: (1) hi Ref. |^, further evidence of an 
anomalous proximity effect in the cuprates is presented. 
These results add to those of Ref. [H EM showing that the 
anomalous effects are real. (2) In Ref. I32L l33l the phase 
diagram of YBCO was recently investigated in the pres- 
ence of Ca doping. Among many results, it was shown 
that a glassy state is generated between the AF and SC 
states in Ca-doped YBCO, with a phase diagram very 
similar to that in LSCO and our Fig. EJd). This re- 
sult suggests that Ca-undoped YBCO may have either 
a region of local coexistence of SC and AF or a first- 
order transition separating them (as in Fig.|7[a,b)), and 
only with the help of extra quenched disorder is that a 
glassy state is generated. Then, the generic phase di- 
agram of the cuprates - which usually is considered to 
be that of LSCO - may not be as universal as previ- 
ously believed, as discussed in this publication. Our 
study showing that bilayered systems are more stable 
than single layers with respect to disorder is also compat- 
ible with the experimental results of Ref. I32L l33l namely 
the 1-layer material is more likely to have a glassy state 
between AF and SC than 2- or higher layer materials. 
(3) Our effort has already induced interesting theoret- 
ical workj^ in the context of J-U models. (4) The- 
oretical work(3^ closely related to our proposed glassy 
state in Fig. [7| has addressed inhomogeneous Josephson 
phases near the superconductor-insulator transition. (5) 
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Recent neutron and Raman scattering investigations ap- 
plied to La2CuO4.05 has shown the coexistence of SC and 
AF phases in this compound. |36| (6) Finally, our results 
have similarities with those recently discussed in the con- 
text of Bose metals as well. 371 



The study also provided predictions for photoemission 
experiments (to be discussed elsewhere) and a simple ex- 
planation for the T c increase with Ng (another explana- 
tion can be found in Ref. l3^) . Clustered states are crucial 
in manganites and other compounds, j3(| and this anal- 
ysis predicts its potential relevance in HTS materials as 
well. 



FIG. 7: Schematic representation of the phase diagrams that 
our models show in the clean (a,b,c) and dirty (d) limits. The 
theory discussed in this paper shows the possible appearance 
of regions with local coexistence of AF and SC (panel a) , or 
a first-order transition separating AF from SC (panel b) with 
the first-order character of the transition possibly continuing 
in the AF-disordered and SC-disordered transitions, or an in- 
termediate striped regime (panel c). Possibilities (a) and (b) 
have already been discussed in Ref. Htl although here we do 
not invoke a higher symmetry group such as SO (5). The main 
result contain in this figure is the proposed phase diagram in 
the presence of quenched disorder (panel d). Shown are the 
glassy region, proposed to be a mixture of SC and AF clus- 
ters, and the T* where local order starts upon cooling. This 
phase diagram has similarities with those proposed before for 
manganites, 0, 13 and certainly it is in excellent agreement 
with the experimental phase diagram of LSCO 
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